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Abstract 

This is a short account of recent joint work with T. Chmaj and A. Rostworowski on asymptotic 
behavior of hnear and nonhnear waves, as presented at the conference devoted to Myron Mathisson 
held at the Banach Center, Warsaw 2007. 
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I. INTRODUCTION 



The work presented in this talk is part of a long-time project aimed at the detailed quan- 
titative description of the process of relaxation to equilibrium for nonlinear wave equations 
defined on spatially unbounded manifolds. By equilibrium we mean here a stable stationary 
solution, like a soliton, a black hole, or just a flat space. The convergence to these solutions 
occurs through a mechanism of radiating an excess energy to infinity. For a large class of 
physically interesting systems the late stages of this process are universal: for intermediate 
times the convergence has the form of exponentially damped oscillations (called quasinormal 
modes) and asymptotically it has the form of polynomial decay (called a tail). This very 
last stage of the relaxation process, the tail, is the subject of my talk. 

The presentation of this talk at the conference devoted to Myron Mathisson is justified 
by the fact our results shed new light on Huygens' principle, one of the main subjects of 
Mathisson's mathematical interests. Recall that a wave equation is said to satisfy Huygens' 
principle if: (i) the solution at a point P depends only on the initial data at the intersection 
of the past light cone of P with the Cauchy hypersurface or, equivalently, (ii) the solution 
vanishes at all points which cannot be reached from the initial data by a null geodesic (i.e., 
there is no tail). A prototype equation satisfying Huygens' principle is the ordinary wave 
equation in d + 1 dimensional Minkowski spacetime for odd > 3. Actually, according 



to Hadamard's conjecture [l] this is the only (modulo trivial transformations) huygensian 
linear second-order hyperbolic equation of the form 

g^''{x)V^V,<P + A^(x) + B{x)<P = . (1) 



Mathisson proved this conjecture in the case of four dimensional Minkowski spacetime 



Counterexamples to Hadamard's conjecture, which have been found later (see [s] and Roy 
McLenaghan's talk at this conference), do not change the fact that Huygens' property is a 
very rare and unstable phenomenon. Thus, it is natural to ask if there are perturbations of 
the free wave equation which preserve Huygens' property approximately, in the sense that 
the tail which is left behind the wave front is very small. The existence of such special 
perturbations in higher even dimensions is a byproduct of our studies of tails. 



fl. 



2 



II. MODEL AND ASSUMPTIONS 



We consider equations of the form 

n(i) + v{x)(i) + N{(j),V(j),x) = o, n = d^-A, {t,x)eR'^+\ (2) 

for spherically symmetric smooth initial data with compact support. Since we want the 
free part to satisfy Huygens' property, we restrict ourselves to odd spatial dimension d> 3. 
Apart from obvious mathematical motivations, there are at least two physical reasons for 
studying higher dimensions d > 3. First, for linear wave equations higher dimensions are 
equivalent to higher spherical harmonics. This follows from the identity 

which relates the I — radial wave operator in d space dimensions with the radial wave 
operator in three space dimensions for the Ith spherical harmonic with I — {d — 3)/2. 

Second, some geometric wave equations in 3 + 1 dimensions are equivalent to scalar wave 
equations in + 1 dimensions for c? > 3. For example, equivariant wave maps from R^~^^ 
into satisfy the following equation 

which, after substitution ip — rcf), becomes the nonlinear scalar wave equation in 5 + 1 
dimensions 

4 \ 4 

— dr\4>-\ — (p^ + higher order terms = . (5) 

r J 3 



III. TOOLS 

In this section we recall two elementary tools from the theory of linear wave equations. 

The first tool is a formula for the spherically symmetric solution of the free wave equation 

d 3 

□0 = in d + 1 dimensions (hereafter we shall use I — — - — instead of d) : 

This solution, which is a superposition of ingoing and outgoing waves, is parameterized by a 
single function a(r) uniquely determined by initial data (the superscript in round brackets 
denotes the /c-th derivative). 



The second tool is the Duhamel formula for the solution of the inhomogeneous free wave 
equation 00 = F{t, r) with zero data 

t t+r-T 
|t-r-r| 

where is the Legendre polynomial of degree / This expression can be easily obtained 

from the standard Green's function formula by integrating out the angular variables J]. In 
terms of null coordinates u = t — p and v = t + p the Duhamel formula ([7]) takes a more 
convenient form 



<j>{t,r) = ^^^^ J dvj F{u,v){v-uy^'P,Wu, p= '-^^ ^. (8) 

\t-r\ —V 

The formulae ([6]) and ([8]) will be used repeatedly below. 



IV. LINEAR TAILS 

For the clarity of presentation we first consider the linear equation with a potential 

D4> + XV(P = 0, (0(0, r), r)) = {f{r),g{r)) . (9) 

The prefactor A, introduced for convenience, will be assumed small and used as a pertur- 
bation parameter. In order to determine the long-time behavior of 0(t, r) we define the 
perturbation series 

= 0o + A0i + AV2 + ..., (10) 

where 0o satisfies initial data and all higher 0„ have zero data. Substituting this series 
into equation we get the iterative scheme 

□00 = 0, 001 = -V^0o, 002 = -V^0i, etc. (11) 

which can be solved recursively using the formulae ([6]) and ([H]). Assuming that V{r) ~ r~" 
{a > 2) for r oo, we showed in jsl that the leading order asymptotic behavior at timelike 
infinity (fixed r and t ^ oo) is given by 



, , , C(La) 



(12) 



where 

^^^'") = -(27TI)!! (2) ^= j a{u)du. (13) 

—00 

The constant A is the only trace of initial data. The expression (1121) was first derived by 
Ching et al. jf| who used Fourier transform methods. 

We claim that the first iterate provides a good approximation of the entire tail if A is 
sufficiently small, that is 

r) - A0i(t, r) ~ C(A2)t-("+2') . (14) 

This basically follows from the fact that all higher-order iterates 0„(t,r) decay in the same 
manner (or faster) as 0i(t, r). Of course, the main issue is whether the perturbation series is 
convergent; for I = (i.e., in three space dimensions) this was proved in f?]) but for higher I 
the problem is open. Note however that for practical purposes it is sufficient that the series 
is asymptotic to the solution. 

The numerical verification of (fl^ shows perfect agreement with analytic predictions j^. 
We remark that numerical simulations of tails are not quite trivial even in the radial case 
because discretization errors generate artificial tails which might mask the true behavior. 
To eliminate such artifacts one has to use high-order finite difference schemes. In addition, 
quadruple precision is needed to suppress the accumulation of round-off errors during long- 
time simulations. For these reasons the simulations of tails are computationally expensive. 

V. NONLINEAR TAILS 

In this section we consider equation ([2]) without a potential. For simplicity, we take a pure 
power nonlinearity with an integer exponent p > 3 (the generalization to other nonlinearities 
is straightforward) 

D(P-r = 0, (0(0, r), r)) = ief{r),eg{r)) . (15) 

x^:=x-{x-l) (.T-(fc-l)), fc>0, 

x'^ := X ■ {x + 1) (x+(fc-l)), fc>0. 



Here we use the notation: 

x^ := 1, 
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This time the amphtude of initial data e plays the role of a small parameter in the pertur- 
bation series: 

= £00 + £^1 +^^02 + ••• (16) 

Substituting (ITBIl into equation (IT^ we get the iteration scheme 

□00 = 0, 00^ = 0^, etc. (17) 
As above, we get 0o using the formula ([6]) and then evaluate 0p using the Duhamel formula 



(IH]). In the limit of timelike infinity we obtain js] 

C{l,v) 



where 



(18) 



C{l,p) = i-iy-^^^—^[il + l)ip-l)-2f and A= j [a^'\u)r du . (19) 

— oo 

The remarks given above in the linear case apply verbatim to the nonlinear case as well; 
in particular the perturbation series (fT6|) is known to converge for / = [2] and is at least 
asymptotic to the full solution for I > 0. 



VI. COMPETITION BETWEEN LINEAR AND NONLINEAR TAILS 

The most interesting situation occurs when both the potential and the nonlinearity are 
present in equation ([2]). Then, each of these terms produces its own tail: 

linear tail ~ ^-("+20 nonlinear tail ~ t^-^^+^'^v . (20) 

Clearly, the tail with slower decay rate is dominant asymptotically, that is 

0(t,r)~t-'^, 7 = min{a + 2/, (/ + l)p-l}. (21) 

For I = this result (without a coefficient) was first proved by Strauss and Tsutaya 9j. 

To give an example of the competition of linear and nonlinear tails, let us consider the 
Skyrme model. Under a spherical symmetry reduction (corotational ansatz), this model 
reduces to the single nonlinear wave equation for the function F{t, r) 

/sin^ F \ 
dtiwdtF) - driwdrF) + sm{2F) + sm{2F) l^^ + {drF)^~{dtFf\ =0, (22) 



where w = r"^ + 2 sit? F. Regular finite energy solutions of (l22l) must satisfy the boundary 
conditions 0) = and oo) = mvr, where an integer m has the interpretation of the 
topological degree of the solution. For m = 1, equation (1221) has a regular static solution 
S{r) called the skyrmion. This solution is linearly stable and plays the role of a global 
attractor, that is, every solution starting from smooth finite energy initial data of degree one 
remains globally regular for all times and asymptotically converges to S{r). The perturbation 
4>{t,r) = y/w{F{t,r) — S{r))/r'^ satisfies the equation 

(^1 - dl -^d^ (p + V{r)(t) +"^4)^ + higher order terms = , (23) 

where the potential V{r) has no bound states and falls off as for r oo. For this 
equation we have (i = 5(/ = l), a = 6 and p = 3, hence from (!20l) 

linear tail ~ and nonlinear tail ~ . (24) 

Thus, the nonlinear tail is dominant [l^ . This example shows that one has to be cautious in 
drawing conclusions about the asymptotic behavior of solutions of nonlinear wave equations 
on the basis of linear perturbation analysis - even for small amplitude solutions the nonlinear 
effects can be dominant. 



VII. ANOMALOUS TAILS 



An advantage of our approach, in contrast to decay estimates in the form of inequalities, 
is that we control the coefficient of the leading order term of the tail. This allows us to 
identify those exceptional cases in which this coefficient vanishes and the decay is faster. We 
shall refer to such tails as anomalous. 

Let us first consider the linear case with the pure inverse power potential near infinity, 
that is V{r) = Xr~" for r > R. Then, it follows from ( |T3l) that C{l,a) oc ^ — - — ^ = if 
a is an odd integer < 2/ + 1, hence there is no tail in the first order. This means that the 
system 

□00 = 0, D^i = -V^o (25) 

is huygensian. In order to find the tail in this exceptional case we need to solve the second 
iteration equation = —V4>i via the Duhamel formula. After a long calculation (which 
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requires the asymptotic expansion of 0i at null infinity) we get (see [5| for the details) 



A + O 



(26) 



where the coefficient D{l,a) is given by a complicated but explicit expression (see Eq.(23) 
in 0). 

Next, consider the pure power nonlinearity = (j)^. It follows from (JTOl) that C{l,p) oc 
[(/ + l)(p - 1) - 2]^ = if p = 2 and / > 1. Thus, in higher even dimensions the first order 
tail vanishes for the quadratic nonlinearity. This implies that the system 

□</)o = 0, 001 = 02 (27) 

is huygensian. As before, the leading order behavior of the tail can be obtained by solving 
the second order equation 002 = 20o0i via the Duhamel formula. The result (see for 
the details) is 

0(t, r) ^ e'Mt, r) ~ , c = (-l)'^^^^^ j ^ [a^Hri)fdn • (28) 

Note that quadratic nonlinearities occur frequently in nonlinear perturbation theory so 
anomalous tails are in fact quite common. As an example, consider the Yang-Mills field 
in four dimensions with the S0{?>) gauge group, so that the potential Aa{x) is the skew- 
symmetric 3x3 matrix A^^{x). For the spherically-symmetric ansatz 

4'(^) = K^^-^;^') 0(^,0, (29) 

the Yang- Mills equation d^F''^ + = 0, where F^p = d^Ap - dpA^ + [A^.Ap], 

reduces to the scalar semilinear wave equation in 5 -|- 1 dimensions 

{dl - dl - Ur)(l) + 302 + r203 = . (30) 

The quadratic term in (l30l) produces an anomalous tail (128!) which is of the same order, 
0{e'^), as the standard tail ffTSl) produced by the cubic term. Combined together they give 
(see for the derivation) 

+ 00 

0(t,r) ^ £302(t,r) ~ e^cr\ c = j a{u)a' {uf . (31) 
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Finally, we remark that although our analysis of tails was restricted to the flat back- 
ground, many conclusions carry over to more general asymptotically flat spacetimes, in par- 
ticular black hole spacetimes. For example, applying the ideas of section VII one can show 



that the massless scalar fleld propagating outside a 
black hole decays anomalously fast as ~ ^-(3rf-5) 



ligher even-dimensional Schwarzschild 
sl- This suggests that the problem of 



asymptotic stability of the Schwarzschild black hole is easier in higher dimensions. 
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